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SUMMARY 


The flow of fluids with and without heat sources and subject to 
body forces between two plane parallel surfaces which are oriented in 
the direction of the generating body force is analyzed under the condi- 
tion that the temperature vary linearly along.these surfaces. It is 
found that a modified Rayleigh number (product of the reciprocal of the 
ratio of specific heats and the Prandtl number Рг and oe modified 


Grashof number Gra) as well as a parameter ጄል = Pr Gra is of 


с 


significance in this problem; where В 28 the volumetric ВАЗЫ 
coefficient, fy is the negative of the X-component of body force рег 
unit mass, d is the characteristic length, and ер is the specific 
heat at constant pressure. Solutions of this problem are obtained in 
terms of "universal" functions which are tabulated for simple applica- 
tion to specific cases. Representative velocity and temperature dis- 
tributions from which detailed study of the heat transfer is made are 
then computed. When the ratio of CK, (where С is related to the mass 


flow) to the Rayleigh number is of unit order of magnitude, the effects 
of aerodynamic or frictional heating can be appreciable. Asymptotic 
solutions (for large values of the Rayleigh number) which render the 
computations simple are also presented. 


Comparison of the results from the method given herein with those 
obtained elsewhere іп an approximate manner for ә. special case sim- 
lating the natural-convection flow of fluids with heat sources in а com- 
pletely enclosed region shows that the approximate method is suffi- 
ciently accurate for problems in which the modified Rayleigh number is 


less than 10* ç 
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INTRODUCTION 


In recent years the transfer of heat to and from enclosed or par- 
tially enclosed regions by means of natural convection or by a combina- 
tion of natural and forced convection has taken on new significance in 
the fields of aeronautics, atomic power, electronics, and chemical en- 
gineering. Most of the information on these modes of heat transfer under 
such conditions is of a semiempirical or specialized nature; relatively 
little detailed information exists for internal natural-convection flows. 
in reference 1 there appears one of the few attempts to determine theo- 
retically the velocity and temperature distributions in detail and hence 
the heat transfer for an internal flow problem of this kind. In that 
reference a solution was found for the fully developed flow of fluids 
with and without heat sources between two long parallel plates with con- 
stant wall temperatures (where one could be different from the other) 
oriented in the direction of the generating body force. The information 
obtained therein is of practical value in connection with fully developed 
flows subject to body forces where the surfaces are maintained at uniform 
temperatures. Lighthill (ref. 2) employs integral methods to study the 
natural-convection flow in tubes with either one end or both ends closed 
and with constant wall temperatures, and in reference 3 approximate 
superimposed free- and forced~convection flows are obtained for short 
channels апа pipes. 


Ав the next step in the study of natural convection or combined 
natural and forced flows in confined spaces, consideration is here given 
to the configuration of reference 1 with the exception that the thermal 
boundary condition specified is that the surface temperatures vary 
linearly along the plates or surfaces. (One surface, however, may be at 
a different local temperature from the other but the slopes of the tem- 
perature distributions on each surface are taken to be equal.) ‘The 
analogous forced-convection problem is treated in references 4, 5, and 6. 
The present problem similates several important physical occurrences of 
this phenomenon; for example, it could represent the case where the out- 
side of the channel formed by the plates is cooled (or heated) by a 
counterflow. Іп addition, the present problem represents a more general 
case than was considered in reference 1, since here the temperature will 
no longe: be restricted to be a function of the transverse coordinate alone 
and, hence, energy convective as well as mass convective effects will be 


included. 


The solution is obtained in terms of functions which depend on only 
one of the several associated dimensionless parameters, and these func- 
tions are tabulated so that specific cases can be easily computed.  Solu- 
tions for pure natural convection and for superimposed natural and forced 
convection are shown to be essentially identical. Representative velocity 
and temperature distributions are also presented, and the effects of 
frictional or aerodynamic heating on the flow and heat transfer are 
discussed. 
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А special case which simulates a completely enclosed region in which 
there is no net mass flow, the walls are at the same temperature, and 
heat is generated uniformly by heat sources is treated, and detailed 
velocity and temperature profiles are obtained from which the heat trans- 
fer is determined. This special case was treated in an approximate 
manner in reference 7. 


Consideration is also given to the problem of convective inversion, 
that is, to the cases where the modified Grashof number changes sign. 
It is shown that convective inversion due to changes in the body force 
direction, to changes in the sign of the volumetric expansion coeffi- 
cient, or to changes in the sign of the axial temperature gradient alters 
the character of the problem, because, it is believed, under these con- 
ditions the flow becomes unstable because of heating from below (see 
pp. 104 to 107, ref. 8). 


ANALYSIS 
Formulation of the Problem 


The study to be made here is that of the laminar fully developed 
flow of fluids with and without heat sources and subject to a body force 
between two plane parallel surfaces open at both ends and oriented in 
the direction of the generating body force (see fig. 1). It is further 
specified that there shall be linear (with equal slopes) temperature 
variations along the walls but that the walls need not necessarily be at 
the same temperature. The flow is assumed to be parallel to the axis of 
the channel (that is, the only nonvanishing velocity component is the one 
in the longitudinal direction) and in addition it is assumed that the 
physical properties (for example, Cp and ы) of the fluids are constants 
and that the essential influence of the density changes on the flow is 
taken into account by the introduction of the volumetric expansion co- 
efficient in the body force term (that is, the other influences of vari- 
able density and the variation of the expansion coefficient with tempera- 
ture are negligible). Discussions of the justification of the assump- 
tions can be found in references 1, 9, and 10. 


Under the conditions stated, the basic equations with body forces 
included expressing the conservation of mass, momentum, and energy (see 
ref. 1) become, respectively, 


90 о (1) 


xia + Se) (2) 


መ 
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oF = 0 (3) 


6:2. 9” Po әт” SUY а 
E E Sl ) “kc (4) 


(See appendix A for a complete list of the symbols used herein.) With 
cognizance taken of equations (1) and (3) the above system reduces to 


ай) 3 ар 

— = т (ве + 55 (5) 
5፲* gape РС ame a Q 
c^ mit ЫЕ - E “Ж (6) 


where the velocity is a function of the transverse coordinate Y only 
and the pressure is a function of only the longitudinal coordinate. 


The body force term in equation (5) can be written as a buoyancy 
term by introducing the volumetric expansion coefficient В in а manner 
similar to that described in appendix B of reference 1. Equation (5) 
then becomes 


ppt 
4,5% e nee - (Ен) 7 


where the subscript wo refers to the surface at Y = О (вее fig. 1). 


The boundary conditions associated with this problem are that the 
velocity at the walls must vanish (the no-slip condition for viscous 
fluids) and that the temperature mist vary linearly along the walls. 
(Note that the latter condition implies that the temperature gradients 
along the walls and hence the axial heat flux along the walls must be 
constant.) In order to satisfy the temperature conditions and equation 
(7), the temperature must be of the form 


T*(X,Y) = AX + T(Y) (8) 
Mathematically, the boundary conditions are formulated as 
U(0) መ U(à) = O (9) 


T"(x,O0) = АХ + ፲(0) (10) 
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т*(х,а) = AX + т(а) (11) 


Substituting equation (8) into equations (7) апа (6), respectively, 
ylelds 


ай; ፀዞ፤ህ 1 (22 

— + — 6 =  [= f 1 
ала 

2 pe A e 

Ө UY p (ӘУ Q _ 

са t +Ë (F) +ተጄ=ዐ (35) 
where 


9 «Т(Ү) - T(O) 


Since the left side of equation (12) is & function of Y alone 
and the right side is a function of X alone, it is clear that each 
side must be equal to a constent. Thus, equation (12) can be written as 


EN pem Ө = С (14a) 
оғ 
1 (ав = 
m (% ü мо) 2d (14) 


The temperature now appears in equations (13) and (14) as Ө which 
is independent of the longitudinal coordinate X, and hence the thermal 
boundary conditions can be written as 


8 (0) 
ө(а) 


О 
т(а) - ቿ(0) = Oy (15) 


To nondimensionalize equations (14) and (13), let 
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У = yd | (16) 


Bf.d 
X ) is the new dimensionless parameter presented 


where Кл = Ргбгд I 


in reference 1, but here the Grashof number is modified because it is 
based on Ad rather than on a temperature difference, that is, 


Gra = "E Hence, equations (14), (13), (9), and (15) become 
u" + T = СКА (17) 
T" - Ra u + (u')2 + оК = 0 (18) 
u(0) = ъ(1) = О (19) 
т(0) = 0 (20) 
шы (21) 


where the primes denote differentiation with respect to у, 
Ба = (1/ү) Prürj is-a modified Rayleigh criterion or number (see p. 105, 


ref. 8), а = Qd/kA is the dimensionless heat source parameter, and 
c =G Ред? /c AK, - Eliminating T between equations (17) and (18) 
results in | 


ulYV - (3፤)ቾ + Ва u - aKa = O (22) 


NACA TN 5141 ^ 7 


with the boundary conditions 


u(0) = u(1) = 0 (23) 
u" (0) = СКА (24) 
u"(1) = (C - 6,4/А4) Ka = መመል (25) 


where mz 1 - e, САЯ. The constant C which appears &s & parameter in 


the boundary-value problem described by equations (22) to (25) merely 
specifies the temperature level (see eq. (17)) of the problem. In order 
to define completely the temperatures and velocities, this constant must 
be related in some way to the physics of the problem. From equation (14b) 
it can be seen that С could be determined from the pressure gradient 
along the channel; that 1s, C is essentially connected with the end 
conditions to which the channel is subject. Since the pressure gradient 
may not be known а priori, in the subsequent section dealing with the 
solution of the present problem C will be related to the end conditions 
by the mass flow in the channel, which remains invarient over the entire 
length of the channel. 


Note that solution of the preceding boundary-value problem will 
yield velocity and temperature distributions for both natural-~-convection 
and combined natural- and forced-convection flows. The forced- 
convection pressure gradient merely alters the magnitude of the constant 
C. A discussion of such a superimposed flow problem under special con- 
ditions is given in reference 11. 


Several interesting observations can now be made concerning these 
equations. First, comparison of equation (22) with the corresponding 
equation in reference 1 shows that they are identical except for the 
third term in equation (22), which does not appear at all in the equa- 
tion in reference 1. This term stems from the convection term in the 
energy equation; hence, in the present problem energy convection effects 
will be included. The energy convection term vanished identically in 
reference 1 because of the assumption that the velocity and temperature 
profiles were independent of the axial coordinate. Second, since the 
convection term appears with а coefficient, another dimensionless param- 
eter (Ra, the modified Rayleigh mumber) is associated with this problem 
and its influence on the results mist be studied. Finally, the condi- 
tions of the problem require that the temperature be of the form 
T* = AX + Т(у) (see eq. (8)), and hence the longitudinal heat flux is 
everywhere constant. 
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The boundary-value problem stated in equations (22) to (25) can be 
written in more convenient form by defining 


v= Rau = aKA (26) 
Hence 
УУ + ва у - 2 በ") = 0 (27) 
v(0) = v(1) = - аҚд (28) 
у"(0) = + Ва CK, v"(1) = + m Ra CK, (29) 


Solutions of the Boundary-V&lue Problem 


Equation (27) is nonlinear (the nonlinear term is due to the fric- 
tional or aerodynamic heating) and therefore, as in reference 1, a method 
Of successive approximations will be сағылды to find its solution. То 
this end, equation (27) is written 2 


тп + Ra Vy - = (v. br a (30) 


n 


where n 40,1 denotes the particular term in the approximation 
у = VQ ተ уу and vi, = 0. Let 


የ=፳ኛ 

| = ፪፻ - 1 

እ =-= o (8) 
_ Ral 20КА 

B 64 

в = 55 
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Equations (30), (28), апа (29) become, respectively, 


viY + BV, = z(v : 2 = | (52) 

vo(-1) = 00) = እ (55) 

፣6(-3) = JNR ደርጊን = መፓዛጾ (Зда) 

ул(-1) = (1) = vi(-1) = vi(1) = O náo (34b) 


where the subscripts now denote differentiation with respect to 1. 


Zeroth-order approximation. - in the zeroth-order approximation, 


e 
A the nonlinear term which is associated with the frictional heating does 
o not appear and the problem then consists of solving the equation 
" viv + Ву. = O (55) 
О О 
` subject to the boundary conditions 
Vo(-1) = v9 (2) = እ (36) 
50-1) = “ላጆ vo (1) = [В (27) 


For simplicity of computations the solution can be obtained in terms of 
symmetric and antisymmetric functions of т, depending only on the peram- 
eter R (or Ra), by setting 


Jim - 1) dim + 1) 
= АУСО + 2 VOL + ( 2 У02 (58) 


where the boundary conditions to be satisfied by the various voj (where 
j = 0, l, 2) are given in the following table: 


Subscript | v(-1) | v) v(1)| v"(-1) 
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Тһе functions Vo must each satisfy equation (55), which has the general 
solution for positive values of the Rayleigh number 


Vo = 84 exp(ern) + ao exp(-ern) + ax exp(grn) + 84 exp(-grn) (59) 


where 
ешісі 
Еші-і1 
and 
r = RENE 
Using the boundary conditions to evaluate the constants in equation (59) 


and expressing the solution in terms of real products of circular and 
hyperbolic functions show that 


L 
Van = —— n (cosh r сов r cosh rq сов ry 
00 cosh?r + cos?r - 1 | 
+ sinh r sin г sinh ry sin гї) (40) 
—— ሚሙ ; 
Yo. = 5 (cosh г sin г sinh ry cos гт 
сов г - cosh г 
- sinh г сов г cosh гт sin rn) (41) 
l 


Van = (cosh г сов r sinh ry sin rn 
Ой ^ cosh^r + сов2р = 1 


= Binh r sin г cosh гт cos rn) (42) 


Note that "00 and үсә are symmetric functions and Voi is an anti- 
symmetric function of т. These "universal" functions "04 are given 


for various v&lues of R& in table i. 
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First-order approximation. - In order to include the effects of 
frictional heating to a first order, let n - 1 in equation (32). Thus 


ደኛ + Bv. = È (va)? = (0) (45) 


where, in view of equation (58), f(n) 1в given by 


ያ(ካ) = 2 A?v! Z L XJ(m + 1)у: ү! , (m - 107,2 
Y ар |ፆ Yoo = Yoo "02 4 01. 
22(п + 1)? 5 : (аб - 1) Е А 
Ы £ "62. + ATM - l)voo ч + 2 ‘ou “o2 
(44) 
The boundary conditions on v, are 
v,(-1) = v, (1) = У1(-1) = vi (1) = 0 (45) 


Once again, to obtain "universal-type" functions let 


2 
уі = А 910 + ÀJ(m + 1)У414 + Xm- 1)” у 


12 
2 2 2 አፍ 
መኳ + 1 Ju - 1) 


In view of the form of equations (43) to (46), the functions Vig» Уту, 
Vig» and уу are symmetric; and v4, and уу are antisymmetric 
functions. Each of the functions үү, (k 20, 1, .. . 5) will bea 


solution of equation (43), but where the nonhomogeneous term is only 
the related part of equation (44) (for example, the first term on the 
right-hand side of equation (44) is associated with the vig solution 


and so forth). 


A particular solution of equation (43) can be constructed from 


(v, = |. ደ(፻) G(n - £)at 
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where С is the Green's function which satisfies the homogeneous part of 
equation (45) and the conditions G(O) = G'(0) = G'(O) = O, G'''(O) «1. 
These arbitrary conditions were applied to yield а simple form for the 
Green's funetion. Hence 


n 


а ЕРГЕ 2% Ге sinh g r(n - Ё) - g sinh e ፓ(ካ - 2) ያ ቼ )6፻ 
О 


(47) 


Note that if the f, is а symmetric function, (vi p 46 also symmetric; 


and if fx is antisymmetric, so is (уік)ір. From the boundary conditions, 
equation (45), it can be seen that vj, must be either symmetric or anti- 
symmetric. Therefore, if fy is symmetric, the complementary solution 
used with the particular solution must be symmetric; and if fk is anti- 


symmetric, the complementary solution must also be antisymmetric. These 
complementary solutions аге, respectively, 


_ p cosh err cosh grn 
(vix) св = Ë cosh er + T cosh gr (48) 
_ ፌ Binh erm sinh gry 
(ues "Ее ተ sinner ini 
Eom 2 [sinker (1 - t) ደ (5)65 (50) 
ሀ 
L 
Р = ы [sinh g >G - £)| т (фа (51) 
O 


The уур solutions сап be written explicitly in real form (by proper 


combinations of particular and complementary solutions) and are given in 
equations (B1) to (B6) in appendix B. Values of v for several R 


lk 
(or Ra) are presented in table II. 


Velocity and Temperature Distributions 


Now that the various VQ, are known explicitly, they can be inserted 
into equation (46) to yield v4. The sum of Хо (as given by eq. (38)) 
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and ገ then forms the solution which includes frictional heating effects 


to a first order. Ву means of the various transformations made, the solu- 
tions of the given boundary-value problem (eqs. (22) to (25)) or the 
dimensionless velocity distributions to the zeroth and first epproximation 
are, respectively, 


K 1/2 1/2 
А | - 1) R& C (m + 1) Ra C 


к,2 
1 А 2 1/2 
Uy = Ба a T vi) + š = ug + 84፪ኗ | Vig aka / C(m + 1) 711 


2. 2 2 
(m - 1) “вас (m + 1) Вас 1/2 
ተ 5Ё Т Vie + " Vis - (m - 1) ова2/ су 


14 


2 
m" = 1)Ваб 
+ 5 | (55) 


Where the 70; (3 20, 1, 2) аге given in equations (40) to (42) and the 
ук (E = O, 1, 2, 3, 4, 5) are given in appendix B. 


ln principle, higher-order approximations could be obtained by con- 
tinuing the procedure described. However, the results become very un- 
wieldy. Therefore, beyond the range of applicability of the zeroth- and 
first-order approximations (that is, in the range of large frictional 
heating effects), the complete boundary-value problem should be solved 
numerically; some discussion of these numerical results relative to the 
zeroth- and first-order approximations will be presented subsequently. 


To determine the temperature distributions, recall that they are 
related to the velocity distributions by equation (17) and that 


Ra 
п ጠጩ 
Yoo = Yoo and Yoo = | Yoo so that 


1 (m + 1) 1/2 
То = СКА - up EX - UE очор t 2(m - 1)Cvoi t 2 На” voo 


(54) 
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and 


„ (m - 1) вас” “Вас” i. | (55) 


where the primes on the u-functions denote differentiations with respect 
to y and on the v-functions denote differentiations with respect to т, 
and the explicit forms of the second derivatives appearing in equations 
(54) and (55) are presented in equations (B7) to (В15) of appendix B. 


Thus, solutions of the original boundary-value problem in terms of 
u and Т are known to zeroth- and first-order approximations, that is, 
neglecting frictional heating and including its effects to a first order, 
respectively. These solutions are, however, given in terms of the param- 
eter C (recall that m is also а function of С through eg. (25)). 
Therefore, to relate С to the physical problem the dimensionless mass 
flow in the channel is defined ав 


M m ugdy (56) 


Neglecting aerodynamic heating, equation (56) becomes (using the zeroth- 
order approximation as given by eq. (52)) 


1 


5 а а, [sin Zr + sinh 2r 
NS uy = gei 41/4 EORR ar ተ сов =. 
0 
Р Ral/4 C(m ተ 1) (біл 2r - sinh 3 (57) 
^/2 cosh 2r + сов 2r 


or 


Pu. 
A cosh 2r + сов 23 B "E a iem 2r + sinh | 
= 2ለ 1/4 3 7 /4\cosh 2r + сов 2 | 
т + ስ in 2r sinh 2r Е, ва? 4. г 


(58) 
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(For large На the expressions in parentheses in equation (58) reduce 

to 1 апа -l, respectively.) Substitution of C ав determined into equa- 
tions (52) to (55) then yields the velocity and temperature distributions 
for any case if M, O47, / Ай, Ва, а, and Кл are known. Note, then, that 


the mass flow is an independent parameter of the problem, and hence could 
be due in part to & forced flow. For & given configuration and fluid the 
last four of these parameters are specified, so thet equation (58) ге- 
lates С to the mass flow M, that is, the temperature level and the 
mass flow through the channel аге related, as is only reasonable in flows 
of this type. 


Asymptotic solutions. - The solutions presented in the previous sec- 
tions are valid (to the proper order of approximation) for all values of 
the parameters of the boundary-value problem. However, from physical 
considerations, it can be seen that in many practical occurrences of the 
phenomenon under consideration, the parameter Ra шау become very large 
(of the order of 104 and higher). It is therefore appropriate to examine 
the asymptotic character of the boundary-value problem. To this end it 
is convenient to write equations (22) to (25) as 


ጁር 
log". (qu + ú - መጭ = (59) 
u(0) = u(1) = O (60) 
u” (0) ==... ኳሽ (3) =m (81) 
where 
u = Ка (52) 


For very large Ra, equation (59) is of the boundary layer type (see ref. 
10). Therefore the velocity and temperature profiles will have very large 
gradients near the walls, and thus the asymptotic solutions will yield the 
velocity and thermal boundary layers with essentially constant conditions 
given by the inviscid solution чт = a/CRa in the center of the channel 


associated with large Ra flows. Hence, expanding the coordinate normal 
to the wall, ав is done in boundary layer theory, requires that 


ዘ 


y = Ral/4ty (65) 


5 = Вагу (64) 
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The asymptotic forms ОҒ equations (59) to (61) are 


ve + - (yy. ሙግ = 0 (65) 

УУУУ ፲ይ “У 1/2 

CRa 
ү(0) = O (66) 
_ 1 for the wall at temperature Ty 

v—(0) =m = 0 (67) 
Jy m for the other wall at temperature Tua 
V(e) = Үсу(о) = 0 (68) 


where the subscripts denote differentiation. Note that the conditions 
expressed by equation (68) replace the boundary conditions аф the second 
wall and require that influences of one wall do not affect the other. 

Thus y can be considered as the coordinate normal to the first (or left- 
hand) wall (that is, the one corresponding to Y = O) and -y will be 

the coordinate normal to the other wall, which is at Y = d. Hence, use 
should be made of the proper part of equation (67) in each solution. 


From equation (65) it can be seen that the frictional heating effects 
will be negligible for large Ка unless CK,/Ra {or CyBfyd/cp) is at 


least of unit order of magnitude. It should perhaps be pointed out here, 
in contradistinction to the qualitative discussions in references 1 and 

12, that the frictional heating is important only if essentially the 

ratio of K (based on any appropriate temperature ) to Ra is of unit 
order of magnitude or larger, as can, in fact, be verified in general. 

Thus it should be noted that the discussions in those references hold 
specifically only if Ra is of unit order or smaller or if Ra does not 
appear as an explicit parameter of the problem (ав in ref. 1, for example). 
For the range of conditions and physical properties of fluids being con- 
sidered, it is unlikely that CK,/Ra will be of unit order for a flow 


generated іп a gravitational field alone with large Ra. Therefore, un- 
less the natural-convection flow is being generated by a body force соп- 
siderably stronger than gravity, the ratio of the volumetric expansion 
coefficient to the specific heat at constant pressure is unusually large, 
or there is considerable forced flow (to increase C), the frictional 
heating effects will not be important for large Ба. (Of course, the 
possibility always exists that some unusual fluid will be employed whose 
physical properties are such that CK, /Ra Will be of such a magnitude 


that the frictional heating effects will be important for large Ra even 
in а gravitational field. Liquids near their critical state may be rep- 
resentative in this respect (see ref. 15)). Furthermore, it can be seen 


СВ-5 


МАСА TN 5141 | 17 


from the same equation that the effects of the heat sources will be im- 


portant only if о/ CRal/? is of unit order of magnitude for large Re. 
This is physically reasonable. 


Since the frictionsl heating effects аге negligible for large Ra, 
unless CK, is very large these effects will be neglected in this sec- 


tion. (А method of successive approximations similar to that described 
in the previous sections could be applied without difficulty to eq. (65) 
if these effects are of consequence.) Therefore, letting 


V=V- о/ сва М e equations (65) to (68) become, respectively, 
V—— + V = 0 (69) 
УУУУ 
V(O) = = 76 (70) 
СБа 
Ууу(0) = m (71) 
(=) = Vo = 0 (72) 


= „75/1 М m sin —— + ————— cos -L— (T3) 


<j 


To find the temperature distribution the second. derivative of equation 
(73) is necessary, and this is given by 


= -y//2 ee т 
Vyy = -e y/N2 —  — sin -ዱ›- - m сов —— (74) 


сна / 6 2 ላሂ 


For large values of the Rayleigh number На, then, the dimensionless 
velocity and temperature distributions are given by ` 


CK 
А fF a 
u = V + ——- 75 
8 1/2 ( 7s (75) 
Ra СБа, ` 


and 


Ta = СКА (1 - V) (76) 
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Solution of Special Case Simulating an Enclosed Channel 


There is considerable interest in the natural-convection flow in a 
completely enclosed rectangular region. Іп reference 7, the natural- 
convection flow of fluids containing heat sources between two parallel 
planes is considered (as is the case here also); but there, in order to 
simulate an enclosed rectangular region, it is further specified that the 
net mass flow be zero and that the walls be at the same temperature. This 
problem is treated in an approximate manner in reference 7 in that the 
velocity distribution was postulated without regard to the equations of 
motion and, hence, it would be desirable for comparison purposes to obtain 
а more exact solution for this special case from the solutions found 
herein. 


For zero net mass flow in the channel the parameter С can be deter- 
mined directly from the other parameters of the problem and equation (58). 
Thus 


Ө 
со = ТЕ A о, cosh 2r + cos 2r ላ/5 sin Zr + sinh er |. 1 
о  2Ad NE: 86/4 Win 2r - sinh ፳2 ва“ \eosh Zr + cos 2r 


(776) 


where the superscript denotes the zero net mass flow and the subscript 
designates the order of the approximation. It is interesting to note 
from this equation that for no internal heat sources (а = О), zero net 
mass flow in the channel can be obtained when C = 0, / 248. However, to 


obtain the solution for the special case simulating flow with heat sources 
in a completely enclosed region with walls of equal temperatures from the 
solutions presented in the previous sections, 9% must be zero in egua- 


tion (77a), and therefore 


e E a (ов. 2r + сов =] 2 (ша 2r + sinh =. a | 
"E ва sin 2r - sinh 2r Ral 4 \cosh 2r + cos 2r 
(тт) 


The velocity and temperature distributions for this special case are then 
obtained by replacing C in equetions (52) апа (54) by co as given by 


equation (77b). Further physical significance of the zero net mass flow 
case сап be inferred from equation (13). Integration of this equation 
over the channel cross section shows that for no net mass flow all the 
heat generated internally (by heat sources and by aerodynamic heating if 
the latter is significant) in a given cross section must be transferred 
to the walls. 
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Computations made for this case including frictional heating yielded ` 
no appreciable deviations from the zeroth-order results. 


Solutions for Case of Convective Inversion 


An interesting aspect associated with the natural-convection 
phenomenon is that the Grashof number can change sign; this implies a 
reversal in the flow direction and is referred to as convective inver- 
sion (see p. 109, ref. 8). The sign of the modified Grashof number in 
the present paper can be changed in one of three ways: (1) by a change 
of the sign of the longitudinal temperature gradient OT'/OX = A, (2) by 
а change in the direction of the generating body force, and (3) by a 
change in the sign of the volumetric expansion coefficient ав occurs 
near the critical state of a liquid (see ref. 15). 


Since the modified Grashof number appears in the parameters Кл 


and Ra connected with the problem considered herein, the effects of 
sign changes of the Grashof number in the solutions should bg studied. ' 
From its definition, KA is proportional to А, fy’, and В” во that 


only the first will alter its sign. Note further, however, ihat the 
modified Rayleigh number Ra is essentially the product of the Prandtl 
number and a modified Grashof number, that is, a Grashof number which 
depends on the product (Ad) of the longitudinal temperature gradient and 
the distance between the plates. Hence, any one of (1), (2), or (3) 
given in the preceding paragraph will lead to a change in the sign of 
Ra. For negative Ra the solutions as given in the previous sections 
do not apply, and hence the foregoing boundary-value problem (eqs. (22) 
to (25)) would have to be solved with negative Ra. These solutions 

can be readily obtained, but it is found that with frictional heating 
neglected these solutions change character with changes in Ra and that 
there exist critical negative values of Ra for which the solutions 
become meaningless. Іп an attempt to explain these unusual results, 
further interpretation of the problem must be made. Reexamination of 
the meaning of negative Ra shows that not only changes in the body 
force direction and sign of the volumetric expansion coefficient but 
also a change in the sign of the longitudinal temperature gradient A 
can lead to negative Ra. If the negative Ra is attributed to the 
last cause, the physical interpretation of the unusual mathematical 
results pointed out becomes clearer, because a negative A implies that 
the fluid is being heated from below and this situation leads to а 
"Rayleigh-type" instability of natural-convection flows due to the 
"piling of heavy fluid on lighter fluid." Analogous interpretations, of 
course, also follow directly for changes in the body force direction and 
in the sign of the volymetric expansion coefficients. Natural-convection 
flows heated from below between horizontal plates have been studied ex- 
perimentally in some detail (see refs. 8, 14, and 15, for example), and 
it was found that the flow does indeed change character (into cellular 
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motion) for certain critical values of the Rayleigh number. Hence, it 
is believed that the critical values of Ra found from the linearized 
analysis (in that the aerodynamic heating was neglected) for negative 
Ra may be analogous to those observed in actual cases. However, since 
this instability leads to these additional complications and should be 
further investigated (perhaps using the true nonlinear eqs.), the case 
of negative Ra will not be treated further herein. It should, how- 
ever, be kept in mind that if in an actual case of the configuration 
considered herein the Ra is negative, the flow and heat transfer will 
not be as predicted in this paper but should be expected to exhibit a 
behavior pertinent to the 'unstable-type" flows. 


RESULTS AND DISCUSSION 
Velocity and Temperature Distributions 
The relations between the actual and dimensionless velocities and 


temperatures as determined from the various transformations in the 
analysis (see eqs. (16)) are 


A 
йр 
U 55. ` (78) 
Ad. 
9 = KA T (79) 


where U and Ө denote the actual and u and Т, the dimensionless 
quantities. For a given fluid, configuration, heat-source intensity, 
and mass flow, the velocity and temperature distributions can be com- 
puted from equations (52) to (55) (for Ва > 0); and for zero net mass 
flow and the walls of the same temperature, by applying equation (77b) 
to equations (52) and (55). These computations will be accurate within 
the limits of the method of solution; that is, for moderate and small 
values of Ra the solutions yield results of reasonable accuracy for 
small СКА, and for large Ra the zeroth-order approximations or, even 


more simply, the asymptotic solutions will give answers valid for all 
СКА. The range of applicability of the various solutions presented 


herein will be discussed more fully subsequently. 


Because the solutions were obtained in the convenient forms (eas. 
(52) and (54), for example) wherefrom the qualitative effects of the 
various parameters associated with the problem can be studied, and since 
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tabular values of the universal functions аге presented to facilitate 
computations for any specific case, no extensive detailed calculations 
will be given covering the entire range of values taken on by the param- 
eters. Representative velocity and temperature profiles were, however, 
calculated for ጀለ = 10, m = -l, l, and 2, а = О, 10, and 100, and 

Ва = 10, 102, 1600, and 104. Ти addition, the parameter С was given 
the value -1 in all the computations except those for the case simu- 
lating flow іп а completely enclosed region. Іп this мау the relative 
influences of the other parameters are just as apparent, but the number 
of computations is greatly reduced. The results of these computations 
аге presented in figures 2 to 9. The contents of each specific figure 
(numbers 2 to 9) are listed in the following table: 


| m [c [m = 20| 102 | 3600 | ፲* | 30 | 100) 160 | 10° 
-1 O | 2(а) | 5(а) | 4(а) | 5(а) б(а) | 7(а) | 8(8)| 9(а) 
10 2(а) | 5(а) | 4(а) | 5(а) 6(а) | 7(а) | 8(а) | 9(а) 
100 2(a)| 3(a) | ፌ() | 5(a) 6(a)| 7(8) | 8(a) | 9(а) 
1 0 5(ъ) | 24(ъ) | 5(ъ) 706) | መርን | 90) 
10 |8; 

100 | 5(b) 
ከ 
10 4(с) | 5(е) 
100 4(с) | *5(e) 


&Ineludes results for u or T. 
PIncludes results for u(2) or т(2). 
CIneludes results for u, or Tg. 


3(b) 5(b) 7(b) 9(b) 


7(b) 90) 
7(6) | 8(e) 
7(5) | 8(<) 


7(6) 


5(ъ) 


For each triplet of parametric values (m, о, Ra) the profiles were com- 
puted with frictional heating neglected (by eqs. (52) and (54) and ае- 
noted by ug and Tg on the figs.), with frictional heating included 


to a first approximation (by eqs. (53) and (55) and denoted by ui 
and Ту), and in several specific cases with frictional heating completely - 


accounted for (by numerical solution of eqs. (22) to (25) using a Card- 
Programmed Electronic Calculator and denoted by u and T). For 
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Ra = 104 the asymptotic solutions (given by eqs. (75) and (76) and 
denoted by ug and Та) are also included for Кд = 10, m= 2, and 

© = О апа 100. The asymptotic solutions were computed from each in- 
dividual wall to the channel center and then faired in to join smoothly. 
From these computations any qualitative trends obtained by examination 
of the solutions can be further substantiated and, in addition, some 
definition of the range in which the frictional or aerodynamic heating 
exerts а large influence can be made. Calculations were also made for 
the special case simulating a completely enclosed region in which there 
is no net mass flow and the walls are at the same temperature by applying 
equation (77) to the appropriate solutions. These curves are given in 
figures 10 and 11. The velocity and temperature profiles (particularly 
for Ва < 104) are qualitatively similar to those determined experi- 
mentally in reference 16. 


Effect of different wall-temperature configurations (m varying 
and heat sources . - From equations (52) to (55) and their 
related universal functions, it can be seen, as expected, that an in- 
crease in the wall temperature parameter m ог an increase in the heat- 
source parameter a results in larger velocities and higher tempera- 
tures. These trends together with that of increasing net mass flow, as 
represented by the area under the u-curves, with m and о can be ob- 
served on figures 2 to 9. It can also be seen from figures 6 to 9 that 
if sufficient heat is generated by the heat sources, the direction of 
heat transfer will be changed. In agreement with the statements made in 
the section dealing with the zero net mass-flow case, note from fig- 
ures 2(a), - 4(а), and 5(а) that if aerodynamic heating is neglected 
for m = -l (since C = -1) and a=0, there is no net mass flow. Іп 
general, the velocity distributions become more symmetrical with the 
larger а (see figs. 2(a), 3(a), 4(a), and 5(а), for example) because 
the heat added uniformly by the heat sources counteracts any asymmetry 
imposed by the wall thermal conditions. 


Effect of the modified Rayleigh number (Ra). - Examination of 
the solutions (eqs. (52) to (55)) shows that the velocities and tem- 
peratures decrease with increasing values of the modified Rayleigh 
number Ra. This trend can also be seen by comparison of corresponding 
curves іп figures 2 to 9, and even by comparing with the curves іп refer- 
ence 1 which are Тос Ва = О. For large Ra it can be seen from fig- 
ures 5(b), 5(с), 9(5), апа 9(c) that the velocity and temperature pro- 
files take on a “boundary-layer form." Asymptotic solutions computed 
for Ва «104, Кл = 10, m= 2, and с = О and 100 are also presented 
on figures 5(с) and 9(c), and these very closely approximate the more 
exact solutions. Hence, for large Ra the asymptotic solutions can be 
employed to yield reasonable results much more simply. For the case 
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where m = -l and а = О, increasing the modified Rayleigh number changes 
the temperature distribution from essentially the conduction profile 

(that is, an almost linear distribution) at Ra = 10; the increased 
effect of the convection is then apparent for the larger values of the 
Rayleigh number. It is interesting to note that for Кл = 10, a= 10, 

and m = -1, 1, and 2, changes in the Rayleigh number сап so affect the 
temperature distributions that the heat flow direction from one or both 
the walls can be altered. (Compare corresponding parts of figs. 6 

to 9.) This point will be more graphically portrayed in the subsequent 
discussion of Nusselt numbers. 


The velocity and temperature distributions (see figs. 10 and 11) 
for the special case considered herein of zero net mass flow and walls 
at the same temperature are not in general apprecisbly &ltered in shape 
by increases in Ra although the velocity peaks vary inversely with the 
Rayleigh number. The shape of the velocity and temperature profiles is 
seen to be qualitatively the same as that assumed in reference 7. 


Effect of frictional heating. - By comparing the profiles presented 
in figures 2 to 9 computed by neglecting frictional heating (denoted by 
the subscript zero) with those computed including the aerodynamic heat- 
ing to a first order (denoted by the subscript unity), the effect of the 
aerodynamic heating on the velocities and temperatures can be studied. 
Numerical solutions obtained of the complete boundary-value problem 
(eqs. (22) to (25)) in which the frictional heating was entirely taken 
into account are also included (with no subscripts) (see preceding table) 
on figures 2(b), 2(с), 4(b), 6(5), 6(с), 8(b), 10, and 11 for comparisons 
with the approximate solutions. 


In accord with the discussion on the asymptotic solutions, it can 
be seen that when Кл 26 small compared with Ra (recall that С = -1 
in these calculations) the aerodynamic heating effects are negligible. 
Since no computations were made herein for K, > 10, the computations 
made for Ra = 1600 and 10* show no deviation between the zeroth and 
first approximations, and these are also coincident with the numerical 
solution (see figs. 4(b) and 8(b)). Hence, in the range K,/Ra<< 1, 
the zeroth-order approximations will yield accurate results; if, in 
addition, Ra>> 1, the asymptotic solutions provide a simple means of 
obtaining the velocity and temperature profiles. Note that the param- 
eter Ka serves merely ав a scale factor in the zeroth-order solutions. 


In the range where K, and the Rayleigh number are of the same 


order of magnitude, the frictional heating affects the results to greater 
or lesser degree depending on the particular amount of heat addition as 
specified essentially by the parameters m and a (figs. 2, 3, 6, 

and 7) and, hence, the first-order approximations should be employed in 
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this range. Comparison of the first-order solutions with the several 
numerical solutions of the entire boundary-value problem shows close 
agreement for the cases computed. Hence, unless the conditions are 
more severe than the most extreme conditions for the range of parameters 
considered herein as represented by KA/R& = 1, ш = 2, and a= 100, 

the first-order solutions will yield results which include the effect 
of aerodynamic heating of reasonable accuracy. In this range where its 
effects are important, this frictional heat, of course, acts just as do 
the heat sources and leads to increased velocities and altered tempera- 
ture profiles and, consequently, different heat-transfer rates (figs. 2, 


5, 6, апа 7). 


For the special case simulating flow in а completely enclosed re- 
gion, the aerodynamic heating did not affect the results appreciably over 
the range of parameters under consideration. ል numerical solution com- 
pletely including the effects of aerodynamic heating was made for 
Ra = 10; it can be seen in figures 10 and 11 that this solution coin- 
cides with the zeroth-order solution, which neglects the effect о? fric- 
tional heating. 


The complete consideration of Prictional heating (as in the numeri- 
cal solutions) for the problem discussed herein, just as for the case 
(essentially Ra = О) reported in reference l, leads to the two results 
(a) that there exists a critical set of conditions beyond which no solu- 
tions exist, and (b) that where solutions exist there are two solutions 
for every set of admissible parametric values. Examples of these second 
ጠር) are presented in figures 2(b), 2(c), 6(b), and 6(с) (denoted 
by ч(2) апа T(2)), and it can be seen that the velocities are more 
than lO times ав large ав the first solutions апа the temperatures are 
much greater than the corresponding first solution temperatures. These 
185% unusual results cannot be predicted from the solutions obtained by 
successive approximations as described herein, but are found from numeri- 
cal solutions obtained by means. of a Card-Programmed Electronic Calcu- 
lator. Аб present the significance of the second solutions is not ex- 
plained, although it is felt that they are intimately connected with the 
unique regenerative action of the frictional heating in natural convec- 
tion. Тһе existence of the critical conditions appears to be similar 
to the thermal choking phenomenon. 


HEAT TRANSFER 
Nusselt Numbers 
The heat-transfer coefficients for the natural-convection phenomenon 


treated here can be expressed in terms of Nusselt numbers. The Nusselt 
number is here defined as - ን жне б 


hà 1 /8፲* 
wu = — = = 
k ል መ. à 


ЧАСА ТМ 5141 25 


where the double subscript signifies that the temperature gradient is 
to be evaluated at either У = 0 or Y = d, depending on the wall under 
ы consideration. 


In terms of the dimensionless quantities, 


-1 (ат 
Nu = K Ma (80) 


The temperature gradient can be found from the zeroth-order solutions, 
апа the Nusselt numbers сап be computed on this basis from 


ከ360 = oe ae ЕС 2r + sin 2r) 
242 Ба, сова r + cos г - 1 


E 
1 
Q 
ርን 1/2 (n- 1/2 
+ C nri) (sinh 2r - sin ЗЕ Š m-1)Ra (sinh 2r + sin 2r) 
2(cos^r - cosh 2.) 
” (81) 
U Е" ——— Қаша 2r + sin 2r) 
242 Ra сова r + сов r- 1 
1/2 -1)Ral/2 
+ СК Ва (sinh 2r - sin 2) + -- (BE l (sinh 2r + sin 2r 
2(cos ғ - совһ т) 
(82) 


where the first subscript denotes that zeroth-order approximation is 
used, &nd the second denotes the wall with which the Nusselt number is 
associated. (Eqs. (81) and (82) are, of course, specifically for Ra 
positive.) Note that these zeroth-order Nusselt numbers are independent 
of Кд. When СКд is of the same order of magnitude ав Ra, the zeroth- 


order approximation has been shown to be inaccurate; therefore improved 
Nusselt numbers can be obtained by using the appropriate T  (first- 
order approximation) solutions and, for the ваше conditions as for 
equations (81) апа (82), the Nusselt numbers сап be computed from 


š Мало = Nuoo + кп «А - coRal/2 (ш + 1X2 + 2 1) Re 


ë 2 2 ደ አፍ 
+ ©-@ + 1) Я, + обы - авы A, - СР Re | (s) 


26 МАСА TN 5141 


K 2 2 
А 1/2 C^(m ~ 1)^R& 
Nu3 ገ = Nuo = ያፋ үл а), = Сода / (m T 1Х12 + 4 05 
642 Ra 


+ С + VA, + оба - ава А, + C=) pen.) (в) 


where the various Q, in equations (83) and (84) are given by equations 
(B43) to (B48) in appendix B. 


Computations of the Nusselt numbers were made over а range of 
values of Ra from equations (81) to (84) for Кл = 10, m= -1, 1, and 
2, and %= O and 10 where again С was taken to be -1, and the re- 
sults are presented in figures 12 and 13 for the wall at Y = O and 
У = ሷን respectively. The figures show that for the wall at Y 20 the 
Nusselt numbers decrease, in general, with increasing Ra, and for the 
wall at Y = d the Nusselt numbers increase with the modified Rayleigh 
number except for the case where m = -1. For Кл = 10 and a= 10, 
the Nusselt number changes sign with increases in Ra because of the 
variation of the temperature profile with Ra, as was previously noted. 


Figures 12 and 13 also demonstrate clearly the effect of the aero- 
dynamic or frictional heating. This effect, in accordance with all that 
preceded, is extremely pronounced for low velues of the Rayleigh number 
(that is, when K/Ra is of unit order of magnitude). 


Plow in an enclosed region. - For the special case simulating flow 
in a completely enclosed region (M = O, m = 1, and a # О), the calcula- 
tions for the temperature profile were extended over a larger range of 
Ка and plotted in figure 14 ав the ratio Ф of the temperature differ- 
ence to that for pure conduction, as was done in reference 7. The con- 
duction temperature difference used in Ф is the channel center-to- 
wall difference subject to uniform heat generation by sources and is 
equal to aX, / 8. The temperature profiles computed in an approximate 
manner in reference 7 are compared with those computed more exactly by 
the method reported in this paper, and it can be seen that the dis- 
crepancy becomes quite apparent for values of Ra = 10* and above. 


The Ny in reference 7 is related to Ra by Ra = = Nr; hence, for 


most liquids Nr and Ra are identical. If ү #1, Y is merely a 
scale factor. 1% can be seen from figure 14 that for Ra < 10*, the 
temperature gradients at the wall are all almost identical, and hence, 
even if the temperature profiles themselves were not identical, the 
heat transfer computed by the two methods would be in reasonable agree- 


НАСА ТМ 3141 27 


ment. The variations of the dimensionless temperature variable 5 
(ratio of center-to-wall temperature difference to that for pure соп- 
duction) ав used in reference 7 with Ба ав given by the two methods 
are &lso compared in figure 15, and hence the quantitative limits of 
the approximate method can be seen. 


CONCLUDING REMARKS 


An analysis was made of the flow subject to body forces between two 
parallel plane surfaces oriented in the direction of the generating body 
force along which the temperature is specified to vary linearly. The 
solutions for natural convection and those for combined natural апа 
forced convection were found to be essentially the ваше. ፲ቴ was found 
that a modified Rayleigh nunber (product of Prandtl and modified Grashof 
numbers) in addition to the parameter Кд was of significance in this 
problem. Solutions for the velocity and temperature distributions are 
given in terms of "universal" tabulated functions. Detailed velocity 
and temperature profiles were computed and it was found that, in general, 
the velocity and temperature differences increase with the wall tempera- 
ture parameter and with additional heat due to heat sources. The veloc- 
ities and temperatures decrease with increasing values of the modified 
Rayleigh number. When the ratio of СКд to the modified Rayleigh num- 
ber is of unit order of magnitude, the frictional or aerodynamic heating 
appreciably affects the velocity and temperature distributions. Asymptotic 
solutions for large Ra are presented which make computations in this 
range relatively simple. For any given set of the parameters, complete 
consideration of frictional heating implied the existence of two flow and 
heat-transfer states and implied that no solution exists beyond certain 
critical values. 


Consideration was given to a special case simulating the natural- 
convection flow of fluids with heat sources in a completely enclosed 
region with the walls at the same temperature. Computations from the 
solutions for this special case demonstrated that an approximate method 
developed in another paper should yield reasonably accurate results as 
long as the modified Rayleigh number is less than 106, Тһе effects of 
aerodynamic heating were found to be negligible for this case. 


Study of the convective inversion aspect of the present problem led 
to the inference that for negative values of the Rayleigh number, addi- 
tional complications arise because of an instability (due to heating from 
below) of the flow which mst be more thoroughly investigated. 


Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 29, 1953 
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APPENDIX A 


SYMBOLS 


The following notation is used in this report: 


longitudinal temperature gradient 

constants in eq. (39); i = 1, 2, 3, 4 

constants defined by eqs. (B14) to (B35); i = 1, 2, 5,...22 
constant in eq. (17) 

constants defined by eqs. (B36); i = 1, 2, 5, ...6 

constant defined by eq. (77b) 

constent in eq. (14b) 

specific heat at constant pressure 

specific heat at constant volume 

constants defined by eqs. (B37); i = l, 2, 3, ...10 
characteristic length (specifically distance between plates) 


constant defined by eq. (50) 


constants defined by eqs. (B38); i= l, 2, 3, ...6 


constant, (i+1) 
constant defined by eq. (51) 


negative of X-component of body force. per unit mass 


Green's function  - ЖЗ 


modified Grashof number, Bf, Ad“ /y° 


constant, (i-1) 


5141 


МАСА TN 5141 


h 


hi 


Nu 


NT 


el 


heat-transfer coefficient 
constants defined by eq. (B39); i= 1, 2, 3, 4 
constant, Ва1/ “сқ,/ 64. 


Ва 
dimensionless parameter, Pr Gr, e 


thermal-conductivity coefficient 


constants defined by eqs. (B40); 1-1,2 
dimensionless mass flow 
constants defined by eqs. (B41); 2 = 1, 2, 5, 4 


constant defined by eq. (25) 

constant defined by eq. (64) 

Nusselt number, Һа/к 

modified Rayleigh criterion as given in ref. 6, үБа 
pressure 


Prandtl number, срЫ/к 


heat due to heat sources 
constants defined by eqs. (B42); i = 1, 2, 3,...8 
constant, Ra/16 


modified Rayleigh number, = Pr Gra 


constant, RY “| ^/2 
temperature 

veloeity 

dimensionless velocity 


dimensionless velocity, u/CK, 
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dimensionless velocity defined by eq. (61) 


ጨ 
dimensionless velocity, V - ECT 2 


dimensionless velocity, v/64 

dimensionless velocity defined by eq. (26) 
longitudinal coordinate 

transverse coordinate 

dimensionless transverse coordinate 
dimensionless transverse coordinate, Ral/4y 


dimensionless heat-source parameter, Qa /kA 


е ) 
coefficient of volumetric expansion, "Ë : | 


ratio of specific heats 
dimensionless coordinste, ву - 1 


temperature difference, T - Tuo 


constant, -aK,/ 64 


absolute viscosity coefficient 
kinematic viscosity coefficient 
dummy variable 

density 

dimensionless temperature difference 


dimensionless temperature difference, BT /oK, 
dimensionless center-to-wall temperature difference, 8(T)y..1/2/oK, 


constants in eqs. (80) and (81); i= 1, 2, 3,...6 
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Subscripts: 

а, asymptotic solution 

са, complementary antisymmetric solution 
св complementary symmetric solution 
I inviscid solution 

n order of approximation 

P particular solution 

WO conditions at у = 0 

"1 conditions at у= 4 

ሀ zeroth-order approximation 

1. first-order approximation 
Superscript: 


(2) 


second flow and heat-transfer state 
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FIRST-ORDER SOLUTIONS 


The explicit forms of the first-order solutions to be used in equa- 
tion (46) аге: 


V19 = —7z[n cosh rr сов rj + Bo sinh гц sin rr + Bs cosh ern cos arn 
R& 
+ Ви sinh 2rn sin ጃካ + Bs(cos 2r« - cosh 2р1) + 15 93 (B1) 
1. У 
улт = Bg cosh г] сов rj + By sinh rn sin rn - B4 cosh ern сов zry 
ЕЁ 


+ Bs sinh 2rn sin 2гп - 15 В (в2) 


1 
12" 572186 cosh ГП сов гп + Bg sinh гт sin гї 
+ Вус cosh 2 rn сов 2 rn + В) sinh 2 rr sin ery 


%В,-(совһ 2rn + сов 2rn) - 15 By ol (B3) 


ጊ. 
Yas [53 cosh rì cos rr + By, sinh ry sin ry 


- В. cosh Zr: сов 2rn -= B, sinh Zr sin 2r) 


+ Bs (cos 2rn - cosh 2х1) - 15 B3] (B4) 


1 
ኛገፈ = -5 5 B. = sinh ry cos тї + B4 


cosh гї sin гт 
в | 


6 


+ Ву; sinh ёгт cos arn ተ Big cosh ern sin ern 


ጉተ 
B4 sin ёгт + B, 


ዕ 


B sinh 2r] | (B5) 
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1 
Vig = EE sinh rj cos гт + Boo cosh гї sin ry 
- ጅ16 sinh ery cos ጃካ + By7 cosh ern sin erq 
+ Boo sin ary - Big sinh zen (B6) 


The constants B,(i = 1, 2, ... 22) appearing in the preceding equations 


are readily computed for a given Ra. Іп the subsequent section of this 
appendix these constants are written explicitly in & form suitable for 
reasonably rapid computation. 


The second derivatives appearing in equations (54) and (55) are 


11 -Ral/ 2 


ул. = ——— RP cosh r sin г cosh гї sin rq 
01 А( сов“ - oL 


CB-5 


+ sinh г сов г sinh ፓኻ cos тт (в7) 


Y 
ፍብ 
-a 


Vio = == [В2 cosh гї сов rj - Ву sinh rr sin үң 
። + “В, cosh Zry сов ጃካ - £Bz sinh 2г sin 2гї 


- 2В5(сов 2r + cosh 2rn)| (B8) 


l 
үз, = Ta (By cosh rj cos rr - Bo sinh ry sin ry 


+ АВ: cosh cry cos ጭካ + АВ, sinh arn sin 2rn| (во) 


1 
У] 2 = ча Вэ cosh гї сов ГТ) - Bg sinh rn sin rq 
+ 4811 cosh ጃካ сов ern - 4816 sinh ጭካ sin arn 


+ 25, „(cosh 2r] - cos 2rn)| (B10) 


1 
Viz = 58: [224 cosh rT сов гї - B43 sinh r sin гт 


- 4B, cosh earn сов ጭካ + АВ: sinh 2гт) sin 211 


- 2В-(сов arn + cosh 2ет)| (в11) 
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via = = -B45 cosh rn sin rj + Big sinh rn cos rq 
- 4Ву7 cosh 2гт sin 2тт + 4В-д sinh ry cos ery 


- ZB4g sin ጃካ + 2820 sinh 221] (в12) 


715 = = -Bo; cosh rr sin rm + Boo sinh ry cos rq 


+ 4B18 cosh 2rn sin 2rn + 4817 sinh 2r cos arn 


- 826 sin 2гт - 2B4g sinh 2rn| (B13) 


Constants for First-Order Approximation 


The constants appearing in equations (В1) to (B13) are written in 
an expeditious form for computing as follows: 


В, = - EC»: +р„)Е,  (D,-SD,)E, + 2(D,:D,)E, + 2(D,+D,)By 
hah 
512 Byhe 
+ (Dz+D4)Es - (21402) Ев + т B - (814) 
Bo = = г [Е + (50- -22)Е2 + 2(D44Dz2)Éz - 2(Dz1D4)E, 
2 
128 (2517 + hp®) 
- (D1+D2)Es - (05+0,)86|- 3 poc (B15) 
b-h 
52 hihe 
Pass MEC (B16) 
2 2 
is (817-829 
В, = 15 — M о (817) 
2 2 
16 (^ + ከ2") 
Bg = - == карс NE (в18) 
В, = m [ n, -D 2)Es + (D; -5р DE + (D, DE, + (D. + D, )E- 


256 2 2 | В 
15M, ( 2 1 ) ( 19) 
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128 እኸ 
phe 
By = г-|(-0550,)8; + (501-0205) - (Dy+Dg)Eg + (Dg:24)Ec] + — m 
(B20) 
4 
B, - i5 (2526 -ZD,)C. + (2D,4D.-D.)C, - 2(D5+Dg)Cz - 2(D,*D4)C, 
448 448 bh, 52(M_4M, ) 
+ (D5-Dg42D4)C5 + (~2D 4D 7-Dg)Cg) + -m s (B21) 
ее 
52 ከከ 
16(Ms-M, ) 
; | 
в. = £ [(-5Фу+р)®, + (53-ጆ94)ጆ> - 2(Dz+D4)Ez - 2(D.+D5)EA 
512 hho 
+ (Dg3D4)E5 - (91-53)ኞ6|- “5: (B26) 
4 
i 


(зни) веј - EA 2 + 290 (327) 
515 = r= [C2147 Pg) + (D5-Dg-2D4)Q2 - (25+06)95 + (Dr:Dg)Q4 


- (Dg sinh 2r)Qs - (Го sinh 2r)Qg - (Dg sin 2г)©у - (Dig sin 2r ) Qe 


7291 5 4 | ( 


22 
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። г. -Dg3Dg42D,)Q, + (2D44D,-Dg)Qo + (Dj4Dg)Qg + (25426) 94 


- (230 sinh 52)95 + (Dg sinh 52)98 - (Dio sin 2r)Q, 


+ (Dg sin 2r )Qg | ray; |125; - ቾ0ኳለ - 6 sinh er + 18 sin 2x] 


(859) 
16h. 

- DE (850) 
с ix (B31) 
8 sinh er ( 
= - ——— 852) 

эз 
.. 8 625 ይሮ (B33) 


= Fy [(-221-27405)92 + (25-06-204)91 + (25*56)9ፈ + (руна) 


+ (Dg sinh 2r)Q7 + (Dio sinh 2r)Qg - (Dg sin 2r)Qs5 


= (Dig sin ar) Q| + төне [204 + 1265 + 18 sin 2r + 6 sinh 2] 


(854) 


ш г-|(9-2,-2,74, + (22,-0,-0,79, - (D4505)Q, + (0.90), 
+ (Dig sinh 2r)Q; - (Dg sinh 2r)Qg - (Dio sin 2r)Qs 


+ (Dg sin 2r)Gg] - там; |125 + 2005 + 6 sin 2r + 18 sinh 2r] 
(B35) 
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- $ sinh Zr sin 2r - $ sinh г sin г 
2 4 
z cosh Zr cos 2r - 2 + = cosh г cos г 
2 2 4 
- Е 668 2r - = cosh Zr + Е cosh r сов r 
(B36) 
4 4 4. 
= = сов Zr + = cosh 2r + = sinh r sin r 
5 5 5 
Š sinh ar sin er + 2 sinh г sin г 
Ё cosh Zr сов or +2 - 22 cosh г cos г 
совһ©г cos?r sinh r sin г 
sinh Š> sin 
cosh?r cos?r 
= sinh^r sin?r cosh r cos r 
= cosh г cos sinh“r 
(B37) 


cosh?r cos г sin@r 


sinh?r sin г сов т 


5 


cosh*r sinr sinh г 


sinh г cos г 


cosh г віп г 
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= cosh Zr сов 2r + 2 - 5 cosh г cos г 
2 8 
= sinh cr sin 2r - = sinh r sin г 
- Š сов 2r + Š cosh 2r - $ sinh х sin r 
(B38) 
4 4 8 
- g cos ег - Е cosh er + = cosh г сов г 
ES 8 
- = Binh 2r sin 2r += sinb г sin r 
2 8 
- g cosh 2r cos Zr ተ 2 - Б cosh r cos г 
cosh r cos r 
sinh r sin r 
(B39) 
cosh r sinh r(1 - 2 coser) 
- gin r cos r(2 cosh?r - 1) 
(cosh“r + cos?r - 1)? 
(B40) 
(cosh?r - cos?r)? (cosh?r + cos?r - 1) 
(cosh?r + cos“ - 1)* 
(cosh?r - cos?r)? 
(cosh?r + ር68ጭ - 1) (cosh?r - cos“) (B41) 
cosh@r sin’r 
віпбү cos“r 
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2 4 
Qi = = cosh zr sin 2r - z cosh r sinr 
2 4 
бо = = sinh er сов 2r - = sinh r cos r 
2 4 
Qz = = sinh Zr cos Zr - = cosh r sinr 
= - 2 cosh ሯ sin 2r + Š sim г сов г 
«7-75 5 
(В42) 
= = 2 sin ar pe sinh r cos r - Š cosh r sin r 
457-75 5 5 
95 = - £ sin Zr + Š cosh г sin г + Š sinh г cos г 
97 = ~ Š sinh er o cosh r sin г - 2 sinh r сов r 
|7765 5 5 
= Š sinh ar = 8 sinh г сов r ~ Ее h r віп г 
96 = = 5 5 °° 


Nusselt Number Constants 


The constants appearing in equations (85) to (84) are given 
explicitly as 


531 = (Bj-B,)sinh г сов г - (Bi4B5)cosh r sin г * 8(B,-Bz)sinh 2r cos 2r 

- 8(Bz+B4)cosh 2r sin 2r - 4Bc(sinh 2r - sin 2r) (ваз) 
922 = (B,-Bg)sinh г cos г - (Bg*B,)cosh г sin г + 8(Bz+B,)sinh 2r cos 2r 

4 8(В,-В.)совҺ 2r sin 2r (B44) 
92, = (Bg-Bg)sinh г сов г - (Вд+Вд)совһ г sin r + 8(B, 1 -B, g)sinh 2r cos 2r 


- (В, +В; )совћ 2х sin 2r + 4В.-(віпһ 2r + sin 2r) (B45) 
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92, = (B44-B42)siuh r сов r = (В)5%В)д)совһ r віп г 


+ 8(Bz-B,)sinh 2г соз 2r + 8(Bz+B,)cosh 2r sin 2r 
- 4B-(sinh 2r - sin 2r) (846) 
(2 5 = (В- 6-ጅ15)6085 r сов r = (815+8 6) sinh г віп т 


+ 8(В-д-Ву+т)совһ Zr cos-2r - 8(B,7+B,,)sinh Zr sin 2r 
- &B4g cos 2r + 4520 cosh 2r . (B47) 


E (Boo-Bo4)cosh r cos r - (ጅ21ቶ5227)5435ከ г sin r 


+ 


8(В1+т+В1д)совһ 2r cos 2r - B(B42-Big)sinh Zr sin 2r 


- 4В20 cos 2r - 4516 cosh 2r | (B48 ) 
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Figure 1. - Schematic sketch of configuration considered. 
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Figure 2. - Dimensionless velocity distributions for various 
heat-source parameters with Ва = 10, к, = 10, C= -1. 
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Figure 2. - Continued. Dimensionless velocity distributions 
for various heat-source parameters with Баж 10, ጅለ = 10, 
С = -1. 
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Figure 5. - Dimensionless velocity | distr tutions for various 
heat-source parameters with Ra = 10°, Кл «10, C = -1. 
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Figure 5. - Continued. Dimensionless velocity dis}ributions 
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C = -1. 


МАСА ТМ 5141 


_ ү 
— dae 
| ИРМА || mA 
Ри ML = 

COAT RNC ET 
_ АРАСЫНА LLL 
Tia A | AAT 
_ NAH 
[LU Му 
LLL LLL NA LI 
НЕСЕТ 


10 


со 


~] 


$, 


н 


Рт 
М 
ИШЕ 


А 
(c) m 


Figure 3. - Conclu de. Dinen в1оп1ев velocity dis tribu tion 
= im в heat-source parameter ае = 102, 10, 


МАСА TN 5141 


Ad 
р 


КА 
с 


Dimensionless velecity, u መ ሀ 


Heat-source 
parameter ን 


(8) ፲6 = -1. 


Figure 4. = Dimensionless velocity distributions for various 
heat-source parameters with Ra = 1600, ጀለ = 10, C = -l. 
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Figure £. - Continued. Dimensionless velocity distributions 
for various heat-source parameters with Ra = 1600, K, = 10, 
C = -1. 
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Figure 4. - Concluded. Dimensionless velocity distributions 
for various heat-source parameters with Ва = 1600, K, = 10, 
С = -l. 
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Figure 5. - Dimensionless velocity distributions for various 
heat-source parameters with Ra = 10 ‚ КА = 10, С = -1. 
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Figure 6. - Dimensionless temperature distributions for veri- 
ous heat-source parameters with На = 10, Ky = 10, C = -l. 
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Figure 6 . - Continued. Dimensionless temperature distribu- 
tions for various heat-source parameters with Ва = 10, 


Ky = 10, C = -l. 
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Figure 7. - Dimensionless temperature distributions for vari- 
ous heat-sourre parameters with Ra = 10", Ky = 10, C = -l. 
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Figure 7. - Continued. Dimensionless temperature distribu- 
tions for various heat-source parameters with Ra = 10", 


Кл = 10, С = -1. 
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Figure 7. - Concluded. Dimensionless temperature distripu- 
tions for various heat-source parameters with Ка = 10’, 
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Dimensionless temperature difference, T = 


Figure 8. - Dimensionless temperature distributions for vari- 
ous heat-source parameters with Ка = 1600, Кл = 10, C = -1. 
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Figure 8. - Continued. Dimensionless temperature distribu- 
tions for various heat-source parameters with Ra = 1600, 


K. = 10, С = -1. 


5141 


CB-9 


МАСА TN 5151 65 


Heat-source 
parameter, 


ACR 

ГАГЫ C LT 
шшк шшш дшш 
оиго 


syi 


Dimensionless temperature difference, T = K, ጊና 


1.0 


(c) m= 2 


Figure 8. - Concluded. Dimensionless temperature distribu- 
tions for various heat-source parameters with Ra = 1600, 


Кл = 10, = -1. 
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Figure 9. - Dimensionless temperature distributions for vari- | 
ous heat-source parameters with Ва = 10^, K, = 10, C - -l. 
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(b) m = 1. 


Figure 9. - Continued. Dimensionless temperature distribu- 
tions for various heat-source parameters with Ra = 10°, 


K = 10, С = -l.a 
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Figure 9. - Concluded. Dimensionless temperature distribu- 
tions for various heat-source parameters with Ra = 10*, 
КА = 10, C = -l. 
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Figure ll. - Dimensionless temperatures for Special case of 
M= O, m = 1, K, = 10, and а = 10. (Note: For Ra = 10, 
the numerical solution « coincides with given curve.) 
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Figure 12. - Nusselt numbers for well 
at Ү-0 and Қ = 10. 
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Figure 15. - Nusselt numbers for wall 
at Y-d and ጀለ = 10. 
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Figure 14. - Ratio of actual temperature difference to that 


for pure conduction for special case where М = O апа m= 1. 
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Dimensionless center-to-wall temperature difference, 


56 = a(t) ተ 1 /8/ ^ 
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Figure 15. - Ratio of center-to-wall temperature difference 


to that for pure conduction for special саге where M = 0 
and m= 1. (Note: Dashed curve is plotted against yRa.) 
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